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Crumpled paper or drapery patterns are everyday examples of how elastic sheets can respond 
to external forcing. In this Letter, we study experimentally a novel sort of forcing. We consider a 
circular flexible plate clamped at its center and subject to a uniform flow normal to its initial surface. 
As the flow velocity is gradually increased, the plate exhibits a rich variety of bending deformations: 
from a cylindrical taco-like shape, to isometric developable cones with azimuthal periodicity two 
or three, to eventually a rolled-up period-three cone. We show that this sequence of flow-induced 
deformations can be qualitatively predicted by a linear analysis based on the balance between elastic 
energy and pressure force work. 



Bending and stretching are the two deformation modes 
of thin elastic plates. For plates of thickness h and typical 
length the ratio between bending and stretching ener- 
gies scales as (/i/i?)^. It results that thin plates, for which 
h <^ will favor bending to stretching if permitted by 
the boundary conditions. For instance, a suspended piece 
of fabric generally exhibits draperies that have a conical 
shape. These particular deformations are almost every- 
where isometric, i.e. the surface is not stretched nor com- 
pressed compared to its initial state. In this example, fold 
number and size are determined by the competition be- 
tween gravity and elasticity [Ij and stretching is focussed 
in a small area in the vicinity of the cone tip. Such coni- 
cal singularities are referred to as d-cone, for "developable 
cone". They are one of the two kinds of elementary singu- 
larities in thin plate deformations, the other being ridges, 
in which stretching is focussed along lines. These two 
elementary singularities can be visualized by unfolding 
a sheet of crumpled paper: they appear respectively as 
crescent- or line-like permanent marks and result from 
plastic deformations of the paper [2j. 

Developable cones have been the subject of specific 
theoretical analyses |3j, as well as experiments consist- 
ing generally in pushing a flat elastic plate into a ring 
[4l|5]. In addition to crumpled paper and gravity- induced 
draperies, d-cones have also been reported in delamina- 
tion processes tissue growth [7], and in suspended 
thin layers of viscous fluid [8j. In this Letter, we propose 
a novel experiment in which d-cones are induced by flow- 
induced loads on a thin elastic plate. In this experimen- 
tal setup, a rich variety of deformations can be observed: 
cylindrical taco-like shapes, d-cones with two- or three- 
fold azimuthal symmetries, and a three- fold cone that is 
itself folded, thereby breaking the azimuthal periodicity. 

The deformation of thin flexible structures due to flow- 
induced loads has been investigated in the past as a pos- 
sible strategy to reduce the drag, comparatively to rigid 
structures [9UTT]. This phenomenon has been particu- 
larly observed for plants that adopt more streamlined 
shapes to withstand high winds or water currents [12 , 13j. 
For instance, tulip trees have broad leaves that can recon- 
figure into cones, thus allowing them to reduce their drag 



and resist breakage This particular reconfiguration 
has motivated analytical [15j as well as experimental [TB] 
studies that have considered circular flexible plates cut 
along one radius, clamped at their center, and placed in 
a uniform flow. These plates have been shown to roll up 
into circular cones becoming sharper as the flow velocity 
is increased. In addition, the drag experienced by these 
plates was shown to be smaller than for rigid bodies. 

In the present study, experiments are performed with 
circular silicon plates of thickness h = 5 x 10~^ m, 
bending rigidity B = 180 x 10~^Nm, density pp = 
1200 kg m~^, and variable radius R. The plate is attached 
at its center to an upstream elbow and held in a horizon- 
tal free-surface water channel of 0.38 m- width and 0.45 m- 
high test section (Fig. [T]). The plate is placed in a steady 
uniform flow of velocity U and density p = 1000 kg m~^, 
perpendicular to the plane that initially contains the 
plate. 

In these experiments, we have varied two control pa- 
rameters: the plate radius R in the interval 0.03 < R < 
0.1m, and the flow velocity U up to lms~^. A typical 
experiment consists, for a given plate of radius R, of grad- 
ually increasing the flow velocity starting from U = 0. As 
the plate deforms due to flow-induced loads, it is imaged 
through the transparent channel walls, both along the y- 
and z-axes (Fig. [T]). Silicon has been chosen as the plate 
material because it allows us to have a flat plate in the 




Figure 1. Experimental setup. A circular silicon plate is 
clamped onto an elbow and placed in a uniform water flow. 
The function (j){0) describe the d-cone deformation of the plate 
as flow velocity U is varied. 



2 



(a) 

ModeC 






Figure 2. Draping modes of a circular plate of radius R = 0.05 m. When the flow velocity is gradually increased, different 
draping patterns become visible starting from "mode C", a cylindrical mode (a), to "mode 2F", a two-fold conical mode (b), 
"mode 3F", a three- fold mode (c), and eventually "mode 3F*", a bent three- fold mode (d). 



absence of flow thanks to the plate rigidity and to the 
silicon density relative to water that is close to 1, such 
that buoyancy is negligible. 

For each plate radius studied, the sequence of defor- 
mation modes is identical. At small flow velocities, the 
plate bends around one of its diameters forming a cylin- 
drical surface with generatrices perpendicular to the in- 
coming flow (Fig. [2^). This mode has been called "mode 
C", for cylindrical. This sort of bending deformation 
has been extensively studied in the literature, both ex- 
perimentally and theoretically. Studies have considered 
three-dimensional geometries such as rectangular plates 
in uniform flow [llj or two-dimensional geometries in 
which a flexible filament was placed into a flowing soap 
film [IIIO]. 

For larger flow velocities, different draping modes ap- 
pear. These draping patterns are shown in Fig. [2]3-d. 
They have approximately a conical shape that can be 
analytically described by the angle (p{0) (Fig. [T]). As the 
flow velocity is increased, a sequence of three different 
conical modes can be observed. Beyond a first thresh- 
old, the plate reconfigures into a cone with n = 2 folds 
(Fig. ^p)' After a second threshold, the plate exhibits 
n = 3 folds (Fig. Both of these deformation are az- 
imuthally periodic such that the function (j){6) has a pe- 
riod 27r/n. This periodicity is lost beyond the last thresh- 
old. The last deformation mode still presents three- 
folds but these folds are themselves "rolled-up" or pushed 
against each other, breaking the periodicity (Fig. ^). 



This last mode persists up to the highest flow veloc- 
ity explored during present study (i. e. U = lms~^). 
These conical modes, which appear in sequencial order 
after the "mode C", are called mode 2F, mode 3F, and 
mode 3F* respectively. It should be pointed out that, 
although the elbow maintaining the plate in the flow in- 
duces a breaking of the azimuthal symmetry in the prob- 
lem, the azimuthal phase of the different modes presented 
in Fig. [2^-d is arbitrary and has been observed to change 
when repeating the experiments. 

To assess the influence of the plate radius on the ve- 
locity thresholds, different experiments have been con- 
ducted with different radii: R = 0.03, 0.05, 0.075 and 
0.1m. Qualitatively, the larger is the plate radius, the 
smaller are the thresholds. The domains of observation 
of the four deformation modes in the plane U-R are re- 
ported in a phase diagram (Fig.[3|. This diag ram corre- 
sponds to experiments performed when the velocity U is 
gradually increased and only reports modes that develop 
spontaneously. Note that, when the flow velocity is de- 
creased, the reverse sequence of modes can be observed 
(i.e. modes 3F*, 3F, 2F, and C), but the velocity thresh- 
olds are systematically lower than for increasing velocity, 
revealing hysteresis loops. These hysteresis loops can be 
as large as 30% of the velocity thresholds and are likely 
due to energetic barriers between modes. 

The draping patterns observed in our experiments can 
be considered as d-cones that result from bending de- 
formations almost everywhere except in the vicinity of 
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Figure 3. Phase diagram in the R-U plane. Each symbol 
correspond to a different deformation mode as pictured. The 
solid lines correspond to iso-values of the Cauchy number as 
defined in Eq. (|3|. 



the cone tip where stretching is localized. To evidence 
this stretching zone, we conducted complementary ex- 
periments with plates made of plastic (PVC). In this 
case, when the plate has been deformed into a mode 
2F, stretch focussing is revealed by two crescent-shape 
permanent marks near the tip. These marks are due to 
plastic deformations, when the elastic limit of the ma- 
terial is exceeded. The same crescent-shape singularities 
have been reported in d-cones obtained by pushing a flex- 
ible circular plate into a ring of smaller diameter [4, 5]. 
Note that, when these plastic (irreversible) deformations 
occur, the plate remains locked in the mode 2F, even for 
large flow velocity. 

To gain better insight into the sequence of draping 
modes observed during experiments, we now consider 
the plate energy associated with the flow-induced con- 
ical deformations. For this purpose, we consider a thin 
circular sheet of bending rigidity 5, radius area S 
and thickness h that initially lies in the plane z = 0. 
When subject to a flow of density p and velocity U in 
the z-direction, this sheet deforms into a conical surface, 
such that the streamwise plate deflection can be written 
C(r, ^) = r/(6>), where (r, are the cylindrical coor- 
dinates with the origin O at the cone tip (Fig. [l]). The 
function f{6) describes the conical modes and satisfles 
/(^)=tan(/)(^). 

In this problem, a Reynolds number. Re = UR/v^ can 
be deflned, which compares inertia to viscous forces, with 
V the kinematic viscosity of the fluid. This Reynolds 
number being of the order of 10^, viscous effects can be 
neglected here. Since buoyancy can also be neglected be- 
cause the plate material has almost the same density as 
the ffuid one, the total energy of the plate £ is the dif- 
ference between its elastic energy and the potential en- 
ergy corresponding to the work done by the ffow pressure 



forces. In the asymptotic limit of thin plates {h/R <^ 1), 
we will assume that the plate undergoes pure bending 
almost everywhere such that elastic energy is dominated 
by bending energy: 



where AC = (/ + f") /r is the local curvature of the plate 
and primes denote derivatives with respect to 0. The cut- 
off radius Rc in Eq. ([T]) comes from the logarithmically 
divergent integral as r goes to 0. Physically this means 
that pure bending deformations cannot persist up to the 
cone tip (r = 0). Around the tip, in a region of typical 
size Rc <^ i?, stretching cannot be neglected anymore. 
Yet, the precise value of Rc is not crucial since the bend- 
ing energy varies only logarithmically with Rc. 

The work done by the ffuid pressure when the plate 
deforms into a cone can be evaluated as 
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where p is the local pressure jump across the plate. With 
the assumptions made above, a dimensionless number 
can be deffned, which compares the typical work done 
by pressure forces (pU'^R^) to the typical bending en- 
ergy {B\n{R/ Rc)). This elasto-hydrodynamic number is 
a Cauchy number and can be written as 



pU^R^ 



B\n{R/Rc)' 



(3) 



This Cauchy number is the only dimensionless number of 
the problem. Iso-values of Cy are reported on the phase 
diagram in Fig. [3] and it shows that the experimentally 
measured thresholds all correspond to a ffxed value of 
Cy' The transition between modes C and 2F occurs at 
Cy ^ 1.3, while transition between modes 2F and 3F 
appears at Cy ^ 40, and the one between 3F and 3F* at 
Cy - 180. 

To go further, we have calculated the energy associated 
with linear draping modes with azimuthal wavenumbers 
n = 2, 3, and 4 (noted modes 2F, 3F, and 4F by anal- 
ogy with experimentally observed modes). In the limit 
of small deffections, the developability condition reads 
lo^ (/^ — /'^) d6> = [2]. The simplest conical surfaces 
with azimuthal wavenumber n that satisfy this condition 
are given by H i| f{0) = /o(l + v'2/(n2 - 1) sinn(9) 
where /o is a small amplitude. To evaluate the pres- 
sure jumps associated to these deformations, we assumed 
that the pressure is constant along a generatrix of the 
conical shape and is equal to the pressure on a circular 
cone of same angle as determined in [TB] using momen- 
tum conservation arguments. It yields a pressure jump, 
p = p/7^(l — sin0). 

To calculate the total energy of cylindrical deforma- 
tions (mode C), we assumed a constant curvature every- 
where on the plate, and a pressure jump, p = pU^^ with 
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Figure 4. Minimal energy as a function of the Cauchy number 
for the modes C, 2F, 3F, and 4F. The calculation for mode C 
is stopped when the radius of curvature reaches 2R/'k. 

Un the component of the flow velocity normal to the plate 
surface. Note that, in order to compare cylindrical mode 
(for which there is no stretching) and conical modes (for 
which stretching is localized in a zone of size Rc)^ we 
used hi{R/Rc) = 4. This value has been considered as a 
fitting parameter and corresponds to Rc ^ 1 mm, which 
is approximately the diameter of the central hole in the 
plate. 

For a given Cauchy number, we can thus calculate the 
deformation amplitude that minimizes the total energy 
for each mode: C, 2F, 3F, and 4F. We can then com- 
pare these minimal energies as the Cauchy number is in- 
creased, or equivalently as the flow velocity is increased 
(Fig. |4|. Although this analysis is limited by the rather 
crude modeling of the pressure and by the assumption 
of linear deflections, the modes that minimize the energy 
follow the same sequence as in the experiments: mode C 
for Cy < 4.6, mode 2F for 4.6 < Cy < 318, and mode 
3F for Cy > 318. Note that the mode with azimuthal 
wavenumber n = 4 never appears to have the lowest en- 
ergy. And indeed, in the experiments, even when we tried 
to force a 4-fold mode, the system always returned to its 
spontaneous two- or three-fold mode. 

The above approach does not allow one to model the 
complex deformations observed for the largest values of 
Cy (mode 3F*). However, for this mode, a qualita- 
tive explanation can be provided. On the one hand, we 
note that the mode 3F* is formed of three folds that 
are grouped whereas they are arranged periodically along 
the azimuthal direction for the mode 3F. Therefore the 
bending energy of the mode 3F* is probably only slightly 
larger than its 3F counterpart. On the other hand, when 
these three folds are grouped, the area swept by the pe- 
ripheral folds is important, and, as a consequence, yields 
a substantial work of the pressure forces. For large flow 



velocities, the importance of pressure work becomes rel- 
atively more important, and thus the mode 3F* becomes 
favored compared to the mode 3F. 

In summary, we have studied how a circular flexible 
sheet deforms when subject to flow-induced loads. A se- 
quence of four deformation modes has been evidenced. 
For low flow velocities, the sheet deforms cylindrically, 
while draping patterns appear for larger flow velocities. 
In this latter case, the sheet deforms into d-cones of dif- 
ferent azimuthal symmetries. As the flow velocity is in- 
creased, a periodic 2-fold mode is first observed, then a 
periodic 3-fold mode, and eventually a non-periodic 3- 
fold mode. The transition between the first three modes 
has been predicted using a linear model based on ener- 
getic arguments and a qualitative explanation has been 
proposed to explain why the azimuthal symmetry is even- 
tually lost for the largest flow velocities. In addition, we 
showed that the only parameter that governs the transi- 
tion between modes is the dimensionless Cauchy number, 
in agreement with experimental observations. 
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